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Abstract 



The effect of the colhnear photon radiation by charged particles is considered in 
the second order of the perturbation theory. Double and single photon radiation is 
evaluated. The corresponding radiation factors are obtained. The QED renormaliza- 
tion group approach is exploited in the next-to-leading order. The results are suited 
to perform a systematic treatment of the second order next-to-leading logarithmic 
radiative corrections to various processes either analytically or numerically. 
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1 Introduction 

The modern high energy physics experiments with advanced techniques and high sta- 
tistical require adequately precise theoretical predictions. Among various effects which 
have to be taken into account, QED radiative corrections (RC) give important contri- 
butions to the predictions. At high energies they are usually computed with help of the 
QED perturbation theory. But direct computations of higher order QED corrections 
to complicated processes can be rather cumbersome. For this reason certain methods 
were developed to evaluate first the numerically most important contributions. In par- 
ticular, besides the expansion in the powers of the fine structure constant a, one can 
use an expansion in powers of the so-called large logarithm, L = ln(M^/m^), where M 
is a large energy scale, and m is a charged particle mass, m -C M. 

In this paper we present the derivation of a particular contribution of QED RC of the 
order O (a^L^'^). It is well known, that the angular distribution of a photon emitted by 
a high-energy particle is peaked in the forward direction. Moreover, it is easy to show 
starting from the matrix element, that a process with emission of a collinear photon 
can be represented in a factorized form (see e.g. Ref. [I]). As usually the factorization 
appear if it is possible to separate the long-distance sub-process of collinear photon 
emission and the main short-distance sub-process. In other words, we assume that the 
experimental conditions of the particle registration allow to neglect the effects of small 
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changes of transverse momenta arising from emission of the photon at the small angle 
with respect to its parent particle: <^ 1. So the cross section (or the decay width) 
of the process with hard collinear photon emission can be represented as a convolution 
of the radiation factor R and the distribution of the radiation-less processes da (in 
example of the 2^2 type): 

dcr[a(pi) + b{p2) c(gi) + d{q2) + j{k ~ (1 - z)pi)] 

= da[aizpi) + b{p2) ^ c(gi) + diq^)] ® R^i'^iz), (1) 
da[a{pi) + b{p2) c{qi) + d{q2) + '^{k ~ (1 - z)qi)] 

= dd[a(pi) + h{p2) ^ c(gi) + d{q2)]R^'^{z), 

where z = E'/ E is the energy fraction of the particle emitted the photon, E and E' are 
the charged particle energy before and after radiation of the photon, respectively. In 
the case of the final state radiation (FSR), we observed the energy of particle c being 
equal to zqi, and we have a direct product of the two factors. In the case of the initial 
state radiation (ISR), we usually compute the kernel cross section in the center-of-mass 
reference frame of particles a{zpi) and 6(^2) and then perform a relativistic boost to 
the laboratory reference frame. 

This paper is organized as follows. In the next Section we re-call the known results 
for the first order collinear radiation factors. In Sect. [3] and in Sect. S] we present 
our results for the second order radiation factors for double and single hard photon 
radiation, respectively. Possible applications of the results are discussed in Conclusions. 
In Appendix we give the explicit formulae for QED splitting functions used in the 
derivation of the factors. 



2 The First Order Approximation 

The derivation of the collinear radiation factors due to an emission of a single hard 
photon in O {a) can be found in Ref. [1]. The factors read 
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The mass of the particle m is assumed to be small compared with the energy, and 
terms suppressed by the factor m? / E'^ are omitted. The photon emission angle with 
respect to its parent particle is restricted by the condition 



< ^0, 



m 



„q2 



(4) 



The energy of the emitted photon is assumed to be above a certain threshold, E^ > AE. 
The parameters "i^o and A either might be related to concrete experimental conditions, 
or serve as auxiliary quantities. In the latter case they should cancel out after summing 
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up the contributions due to emission of the colhnear hard photons with the ones of non- 
coUinear hard photons and of soft photons. These O [a) radiation factors are universal 
and describe colhnear single photon emission for various high-energy processes [1]. 



3 Double Hard Photon Radiation 

In paper [2] the effect of the double hard photon radiation in the Bhabha scattering 
was considered. In particular the effect of the two photon emission inside a colhnear 
cone along the direction of motion of any of the 4 charged particles in this process was 
presented in a form being differential in the energy fraction of both the photons. So 
to get the colhnear radiation factor, we have just to integrate over one of the energy 
fractions keeping their sum fixed. The lower limit of the integral over the photon energy 
fraction is chosen to be equal to the parameter A because both the photons should be 
hard and have therefore energy above AE. In this way for the case of the initial state 
radiation we got 



+ 2/o)f^^(21n(l - 2) - 2 In A - In + Inz 
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+ (1 - z)(2\n{l - - 2 In A - Inz + 3^ + ^-^\n^z^ 



(5) 



where z is, as in Eq. ([2]), the energy fraction of the charged particle after the emission 
of the two photons. 

The corresponding radiation factor for the final state radiation case can be obtained 
from the ISR one by means of the Gribov-Lipatov relation: 
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+ j^(^5\n' z - 2\n z - A\n{l - z) + 41n A + 8 lnz(ln(l - z) -InA)^ 
+ (l-z)f21n(l-z) -21nA-31nz + 3') + + In^ ^1 (6) 



Note that the additional interchanges in the above relation applied for In A and /q 
appear in our case from the crossing relations of the two channels with the given cuts 
on the energies of the soft photon and on the photon emission angle. 



4 Single Hard Photon Radiation 



We have to consider also the process of single hard photon emission accompanied by 
the one-loop virtual correction or by the emission of a soft photon. As concerns soft 
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photon radiation, its contribution does factorize with respect to the coUinear hard 
photon emission: 



dans = Rn ^ Ssda^'\ 5s = ^ , (7) 

where Sg is the one-loop soft photon radiation factor for the given process, computed 
in the standard way [3]. This quantity has an infra-red divergence, which cancels out 
after summation with the virtual loop contribution. And da^^^ is the Born level cross 
section. 

So we would like to get the radiation factor -Rh(s+v)(-2)) where both the one-loop virtual 
correction and soft photon radiation are taken into account. To find this radiation 
factor, we will exploit the known result for the complete second order NLO QED 
corrections provided by the renormalization group approach [1][5|6|7||8] . In analogy to 
QCD we can write the master formula for the corrected cross section e.g. for Bhabha 
scattering in the form (see Ref. [H]): 



1 1 

da = J dz, J dz2Vf:{z,)Vf:{z2) (d(T(°)(-2i, ^2) + dcx^'\z,, Z2) + O {o'L'^)) 




yi y2 



where da^^^ is the O (a) correction to the massless scattering, calculated using the 
MS scheme to subtract the mass singularities. The energy fractions of the incoming 
partons are Zi^2, and Yi^2 are the energy fractions of the outgoing electron and positron. 
Vf^^^"^^^ are the structure (fragmentation) functions of an electron. Here we consider 
only the photonic contributions to the non-singlet part of these functions. The radi- 
ation factors corresponding to the coUinear emission of light pairs were evaluated in 
Ref. [To]. With help of the master formula we can find the most important contribu- 
tions reinforced by the large logarithm L in radiative corrections to a wide class of 
other processes as well. 

We are going to drop the pair contributions, so we need here the pure photonic part 
of the non-singlet structure (fragmentation) functions for the initial (final) state cor- 
rections. These functions describe the probability to find a massless (massive) electron 
with energy fraction z in the given massive (massless) electron. In our case with the 
next-to-leading accuracy we have 

Vf:''^'{z) = 6{l -z) + ^S^)(z,fio,m,) + ^LP^'\z) 

ZTT ZTX 

+ (|^)'Ql2p(°) ® P(°)(^) + LPW ®rf«(^,/io,me) (9) 
+LPg^^>'^^''^{z)^+0(a'L\a'). 

The difference between the functions appear only due to the difference in the next-to- 
order splitting functions P^^'"^), given in the Appendix together with the other relevant 
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functions. The modified minimal subtraction scheme MS is used. We have chosen the 
factorization scale equal to E, and the renormalization scale /io will be taken equal 
to nif.. More details on the application of the approach to calculation of second order 
next-to-leading QED corrections can be found in Refs. [5 ,,8]. 

Let us consider the O (a^L") {n > 0) radiative corrections to a given process, which 
are related to at least one hard photon emission. They can be separated into four parts 
according to their kinematics: 

.(2)NLO_.(2) .(2) .(2) .(2) . . 

"Hard — '^HH(coll) ^ "HH(s-coII) ^ "(S+V)H(n-coll) ^ "(S+V)H{coll) ' \^^) 

where 5hh(co11) gives the contribution of double hard photon emission considered in the 
previous section. The case when one of the photons is emitted at large angle {{}^ > -^o) 

(2) 

and the other one is coUinear is denoted where "s-coU" means a semi-coUinear 

kinematics, see Ref. [2] for details. The term ^(s+v)H(n-coii) corresponds to single hard 
non-collinear {i)^ > ■&()) photon emission accompanied by the O (a) soft and virtual 
photonic corrections. Note that since the non-collinear photon emission doesn't give rise 
to the large logarithm, we can keep in ^(s+v)H(n-coii) ^^^y leading logarithmic terms 
in the sum of soft and virtual corrections. And the last term is the contribution that 
we are looking for: the one due to single hard coUinear photon emission accompanied 
by O (a) soft and virtual corrections. 

From the other hand, the same quantity can be found in the master formula ([8]): 
= ^LPi^' « dag) + I^LPf « da« + ^(LP^^^ + d«) « da^) 



+ (^) (i^p(0)^p{0)+^p(0)^^{l)^^p(l,7)str\ ^^^(0)^ ^^^^ 



2 



2nJ \ 2 



e 



where we leaved out the splitting functions arguments for short. Here daQ^ is the 
contribution of single hard photon emission and da^'* is the soft-virtual contribution 
(in the MS scheme with massless electrons). Lower indexes e and a mean here the 
parts of the corresponding functions related to hard and soft plus virtual radiation, 
respectively. Again we kept in the above equation only the terms reinforced by the 
large logarithm L. 

Comparing the two expression ( fTOl) and ( fTTl) we get 

'^(slv)H(coii) = Rm+v)i^) ® da(^) = £lp(°) ® d^W + |-Lp(r^ ® da« 

+ (—y( — P(0) ® P(0) + LP(°) ® + LP^^^^>A ® da(°) 
V27r/ V 2 /e 



^) Pf^®daW. (12) 
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The MS subtraction leads to the following relations: 



da« = d^l, + da^l - ^{LP^^^ + d^lW\ (13) 
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ZTT 



(14) 



Summing up the parts in (|T2l) proportional to a^^ with help of (fT^ we arrive at 
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(15) 



After substitution (fT3|) and (fT5|) to (fT2l) we get the result 



^H(l+v)(^) ® da(;2) = (— )" LP(°) ® Pf ^ ® da(°) 
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(16) 



Using the tables of convolution integrals [llj we obtain the answer for the ISR factor 



2^^^ (^Li2 (1 - ^) - ln(l - z) In z^ 



-{l + z) z + {1 - z)lnz + z 



^da^'\z), 



(17) 
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(1) 

(S+V) 



•^^Soft + '-1'^' 



Virt 



da(o) 



To get the final state radiation factor we use again the Gribov-Lipatov relation and 
get 



(1) 



jFSR, 



r(0) 



1 + z' 



1-z 

+ (1 + z) In^ ;z - (1 - ;z) In z - 1 



-2Li2 (1 - 2) - 3 In^ z + 2 ln(l - z) In z) 
da(0). 



(18) 



5 Conclusions 



In this way, we received the explicit expressions for the radiation factors, which describe 
hard collinear photon emission in the second order of the perturbation theory within 
the next-to-leading logarithmic approximation. These factors are universal. They can 
be used in analytic and numeric calculations of QED radiative corrections to a wide 
range of processes. In particular, we are going to implement them into the Monte 
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Carlo event generators LABSMC [12], SAMBHA [13], and MCGPJ [H] for several high energy- 
processes. Our results can be exploited also to provide advanced theoretical predictions 
for experimental observables with so-called tagged photons, when hard photons emitted 
at zero (small) angles with respect to colliding charged particles are detected [1^16|17] . 
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Appendix. Explicit Formulae for QED Splitting P\inctions 

The QED splitting functions corresponding to photonic corrections in the leading log- 
arithmic approximation in the first and second orders read 



"l + z^- 


= lim < 


l-z _ 





lim{5(l - z)pf + 6(1 - z - A)p(°)(z)}, 



Pi°^ = 21nA + 



(Al) 



\l + 

+6(l-z-A)2 ^ 



21nA + ^V-4C(2) 



3 \ 1 I ^ 

, 21n(l - z) - \nz + - ) H -\nz- 1 + z 

I - z \ ^ ' 2/2 



where symbols 5 and 6 denote the Dirac (5-function and the step function, respectively. 

The space-like (ISR) and time-like (FSR) next-to-leading terms of the QED splitting 
functions for photonic corrections can be cast as [SfS] 



Pt^>'\z) = 8{l -z)(l- 3C2 + 6(3) + 1±^( -21nzln(l - z) 



+ z + 2U2{l- z] 



1 



1 - z 



{l + z) ln2z + 21nz-22 + 3. 



pa,7)frg(^) = 5(1 _ ^) ('I _ + eCs) + ( 2 In ^ ln(l - z) 



-2\n^z- 2U2 (1 



1 



1 - z 



(l + z) \n^z + 2z\nz-3z + 2. 



(A.2) 



In the next-to- leading calculations we need also the initial condition for the structure 
and fragmentation functions at a certain scale /xq: 



z,Ho,me 



l + z' 
l-z 



In 



l4 

ml 



21n(l -z)-l 



(A.3) 



The dilogarithm and the Riemann zeta-function are defined as usually: 
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Li.(x) = /d,Mi_^^, an) = t^. (A.4) 

n ^ k=l 
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